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We have studied the tunability of the interaction between 
permanent dipoles in Bose-Einstein condensates. Based on 
time-dependent control of the anisotropy of the dipolar in- 
teraction, we show that even the very weak magnetic dipole 
coupling in alkali gases can be used to excite collective modes. 
Furthermore, we discuss how the effective dipolar coupling in 
a Bose-Einstein condensate can be tuned from positive to neg- 
ative values and even switched off completely by fast rotation 
of the orientation of the dipoles. 

PACS number(s): 03.75.Fi, 05.30. Jp 



In the Bose-Einstein condensates (BECs) created thus 
far [Q atoms interact essentially only at very short dis- 
tance. At low temperatures, the interaction between 
these neutral atoms can be typically characterized by 
a single parameter, the s-wave scattering length. The 
magnitude and sign of the s-wave scattering length can 
be changed by tuning the external magnetic fields in the 
vicinity of a so-called Feshbach resonance Q . This tech- 
nique has opened new possibilities for the study and ma- 
nipulation of BECs. In particular, the value of the s-wave 
scattering length can be varied in real time, allowing, 
for example, the observation of collapsing and exploding 
BECs §. 

In recent years an increasing interest is devoted to 
dipole-dipole interactions in ultracold gases . Dipo- 

lar interactions would largely enrich the variety of phe- 
nomena to be observed due to their long-range and 
anisotropic character. If atomic dipole moments are suf- 
ficiently large, the resulting dipole-dipole forces may in- 
fluence, or even completely change the properties of Bosc 
gases [I | 11 1, the conditions for BCS transition in Fermi 
gases [ 7 ] , or the phase diagram for quantum phase tran- 
sitions in ultracold dipolar gases confined in optical lat- 
tices 0. The interplay of short-range scattering and 
long-range interaction may give rise to phenomena like 
ferromagnetic order and spin waves 0. Moreover the 
dipolar particles are considered to be promising candi- 
dates for the implementation of fast and robust quantum- 
computing schemes Q . 

In this Letter, we consider the case of permanent (mag- 
netic) dipoles. For alkali atoms with a magnetic moment 
m of fiB, where /i^ is a Bohr magneton, the long range 
part of the magnetic dipole-dipole interaction is generally 
neglected since it is very small compared to the s-wave 
pseudo-potential. Here, we show that even in alkali BECs 
this dipolar interaction can be made visible if the dipole- 
dipole interaction is modulated by rotating the atomic 



dipoles. For appropriate rotation frequencies, a coupling 
to elementary excitations can be achieved leading to a 
dynamic growth of the amplitude of the excitation. In 
particular, we will show how the dipolar coupling in an 
alkali condensate can be used to excite the quadrupole 
mode and how it is possible to extract quantitative in- 
formation about the ratio of the dipolar coupling and the 
usual s-wave coupling. 

In ensembles of atoms with a larger magnetic moment, 
such as chromium where m = 6/ib, the magnetic dipole- 
dipole interaction energy may become comparable in size 
with the s-wave contribution to the mean-field energy. 
Therefore, even the static properties of a condensate such 
as the aspect ratio are significantly altered. Here, we 
show how the effective strength of the interaction and 
thus its effect on the static properties of a condensate 
can be varied by fast rotation of the atomic dipoles. In 
the case of a dipolar interaction which is stronger than 
the s-wave interaction only certain combinations of trap 
geometry and orientation of the dipoles lead to a stable 
(or metastable) condensate |||]]. By tuning the dipo- 
lar interaction unstable condensates can thus be trans- 
formed into stable ones and viceversa. In this paper, 
we study the interaction between permanent magnetic 
dipoles. However, the methods we propose can also be 
adapted for particles with electric dipole moment such as 
heteronuclear molecules. In this case the dipolar interac- 
tion is expected to dominate over the s-wave interaction 
and thus the tuning of the interaction might be necessary 
to reach stability. 

Manipulating the magnetic dipolar coupling in a 
BEC — The long-range part of the interaction between 
two magnetic dipole moments mi and m 2 of atoms lo- 
cated at ri and r 2 takes the familiar dependence 
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as a function of the interatomic distance r = r 2 — rj., 
where f — r/r and fio is the magnetic permeability of 
the vacuum. We consider the case of a general time- 
dependent homogeneous magnetic field 



B(t) = Be(t) , 



(2) 



slowly-varying with respect to the Larmor frequency 
LdLarmor = mB/h, where m — |mi| = |m 2 | for identi- 
cal atoms. A spin-polarized atomic ensembles will adia- 
batically follow the external magnetic field B(t) and the 
resulting interatomic energy becomes 



U dd (r,t) = - 
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and is thus time- dependent. 

Our analysis is based on the mean-field approach in 
the Thomas-Fermi limit already used in the context of 
dipolar Bose gases fl-|l2]]. Such a description can be ac- 
complished through the hydrodynamic equations for su- 
pcrfluids M: 
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where n is the density, v the velocity field and Sfi the 
difference between the mean-field plus external fields and 
the chemical potential /i given by 



5/^1 = gn + $ 



dd 
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Here M is the atomic mass, MuJqT ■ /2 is an isotropic 
harmonic potential which could be created using optical 
dipole forces, g = Airh 2 a/M and a is the s-wave scattering 
length which we assume to be positive. <&dd(r, t) is the 
mean-field potential generated by the (time dependent) 
dipolar interaction 



$ dd (r,i)= / dr'U dA (r-r',t)n(r') 



(J) 



In the Thomas-Fermi limit it is useful to define a dimen- 
sionless quantity 



£dd 



u m 2 M 
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which is a measure of the strength of the dipole-dipolc 
interaction relative to the s-wave scattering energy. Mod- 
ifications in the ground state condensate density can only 
be expected if £dd is appreciably different from zero. 

For chromium the s-wave scattering length ao is not 
known. If we assume that ao is equal to the sodium 
scattering length (a = 2.8 nm) the dipolar strength pa- 
rameter becomes e£j = 0-29 (m — 6/is). Therefore, 
the dipolar force is expected to lead to visible modifi- 
cation in the density of a chromium BEC. Differently, 
in both rubidium and sodium the magnetic dipolar en- 
ergy is rather small compared to the mean-field s-wave 
scattering energy. The values of £dd are = 0.0064 

and e^l = 0.0035, in the doubly spin-polarized ground- 
state, where m = /is- For heteronuclear molecules with a 
permanent electric dipole moment of one Debye the cor- 
responding dipolar coupling will be increased by approxi- 
matively the inverse square of the fine structure constant. 
This implies that £dd Jlj| is on the order of 10 2 . 

Probing the quadrupole modes in sodium and rubid- 
ium — The time modulation of the dipolar coupling (|^) 
can be used to resonantly excite a quadrupole mode even 
for sodium and rubidium condensates. For simplicity 
we restrict ourselves to the case of a condensate in the 



Thomas-Fermi regime confined in an isotropic harmonic 
trap with frequency ujq (for example ojo/27t = 100 Hz). 
Because £dd is small the dipole-dipole mean-field poten- 
tial <i>dd can be calculated using the symmetric equilib- 



rium density distribution no 
R is its radius 



Muj 2 {R 2 -r 2 )/2g, where 



$dd(r,f) 
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[3 (e(t) • r) 2 - l] r 2 r < R, (9) 



■f^M [3(e(i).f) 2 -l]^ r>R. (10) 
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The potential outside the condensate region (r > R) cor- 
responds exactly to the field generated by N dipoles lo- 
cated in the center of the condensate. Differently, inside 
the condensate (r < R) the mean-field potential has an 
anisotropic but harmonic dependence. 

If the interaction is driven on a resonance frequency 
the evolution of the condensate can be modified even for 
very small dipolar coupling. Consider a magnetic field 
rotating in the plane x — y as B(t) = Be(t), where 



e(t) = cos(i7i)x + sin(£li)y . 



(11) 



The resulting dipolar interatomic energy (|^) and its cor- 
responding contribution to the mean-field potential (||- 
|l(]), will have a 'rotating' part Ej. For small dipolar 
coupling within alkali gases the non oscillating part of 
(^) can be neglected and the time-dependent part will 
resonantly couple to the 'transversal' quadrupole mode 
when u>o/V2 (see Rcf. Jl5[). After linearization of 
the equations of motion (||) and (|^) we obtain 



d 2 Sn 
dt 2 



= V(^V(H 



vO* dd )) 
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where Sn = n — no is the density displacement around the 
equilibrium density no. The time-dependent part of <l>dd 
corresponds in Eq. ( |l2| ) to the rotation of an anisotropic 
harmonic trap with deformation 



J x 
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where ojx and ujy are the maximum and minimum values 
of the rotating trap frequencies. Neglecting the time- 
independent part of (|^) the density will evolve as 



Sn 



3£ d dA^o 
10 5 
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where <f> is the azimuthal angle (r < R). It is noteworthy 
that the rotating part Q of (||) couples onl y t he first 
quadrupole mode in the Thomas-Fermi limit p5| . From 
the study of the evolution of the condensate deformation 



(x 2 + y 2 ) 



r £dd sin(2£!i) 



(15) 



we can extract information about the ratio of the dipolar 
coupling versus the standard mean-field energy. If the 
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FIG. 1. In rubidium the long-range dipole-dipole interac- 
tion does not modify in an appreciable way the ground state 
density. When the magnetic dipoles are driven on resonance 
with a transversal quadrupole mode the anisotropy a (Eq. 
(15)) is growing with time (here we assume a quality factor 
Q = 50). 

quality factor Q = ^/2ujqt of such a resonance, where 
t is the corresponding decay time, is assumed to be 50 
Jil| the amplitude of the oscillation of the condensate 
anisotropy a is expected to be \a\ — 3Qedd/10 ~ 9.5 x 
1CP 2 for rubidium (see Fig. 1). We have performed our 
calculation in the linear regime and for an isotropic trap. 
The non-linearity of the mean-field is expected to give a 
significant correction for larger Q. The extension of the 
calculation to an axially symmetric cigar-shaped trap is 
expected to yield similar results. 

Tunability of the magnetic dipolar interaction: — We 
consider a condensate in an axially symmetric harmonic 
trap, with radial and axial trap frequencies Lu r and ui z . 
The magnetic field 

B(t) = B [cos(<^)z + sin(</j) (cos(fti)x + sin(fit)y)] (16) 

is a combination of a static magnetic field B z directed 
along the z direction and a fast rotating field B p in the 
radial plane. The frequency is chosen such that the atoms 
are not significantly moving during the time while 
the magnetic moments will follow adiabatically the exter- 
nal field B(t); this corresponds to oJLarmor f2 S> utrap, 
where LOtrap is any of the external trap frequencies. In this 
limit we can consider the average of the interaction (^) in 
the period 2 7r/f2, resulting in the cylindrically symmetric 
interatomic potential 

, J , Horn 2 / 3 cos 2 as — l\ , , „ N , 

(tfdd(O) = ( y — J ( 3 cos 6 - !) • ( 17 ) 

The averaged interaction energy ( |l7|) equals (3 cos 2 <p — 
l)/2 times the interatomic energy of dipoles aligned along 
the z-axes. Note that this factor can be changed continu- 
ously from 1 to —1/2, providing the possibility to change 
the dipolar interaction from attractive to repulsive. At 
a particular angle ipM = 54.7° the dipolar interaction 
averages to zero. This angle, the so-called magic-angle, 
is well known in solid-state nuclear magnetic resonance 
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FIG. 2. Tunability of the magnetic dipole interaction. The 
atomic dipoles follow adiabatically the total magnetic field 
( p^[ ) when the field is sufficiently strong. The angle Lp between 
the precessing total magnetic field and the z-axis character- 
izes the interaction, ip = 0: the magnetic dipoles are simply 
polarized in the z direction, (p — tt/2: the magnetic dipoles 
follow the rotating magnetic field and as a result of the fast 
rotation the averaged interaction becomes —1/2 times the in- 
teraction in the polarized case, ip = 54.7° (magic-angle): the 
dipolar interaction averages to zero. 

(NMR) technology The tunability of the dipolar in- 
teraction will allow to observe changes in the condensate 
eccentricity which should be a 10% effect in Cr assum- 
ingacr = aNa or in the collective excitation frequencies 

SH- 

Stabilizing strongly interacting dipolar gases: — Let us 
now apply this tunability to the stability diagram of dipo- 
lar BECs. We consider the case in which the dipolar cou- 
pling is larger than the s-wave pseudopotential (edd !)• 
Here the tunability can stabilize strongly interacting 
dipolar gases like heteronuclear molecules which other- 
wise exhibit instability over a wide range of parame- 
ters. This regime, considered also in Ref. may be 
achieved in a cloud of magnetic dipolar atoms, by making 
use of a Feshbach resonance to reduce the s-wave scat- 
tering length. 

In a spin polarized cloud the stability diagram as 
well as the condensate anisotropy can be studied in the 
Thomas- Fermi limit as function of £dd (see Ref. [^§|). 
When edd > 1 the cloud may be unstable against col- 
lapse, i.e. can be always possible to find a density config- 
uration of arbitrary large negative energy. The stability 
as edd and the trap anisotropy I — w 2 /w r are varied can 
be investigated in analogy to Ref. Q. 

Within the Thomas-Fermi approximation, the stability 
analysis can be made using a scaling variational approach 
|]l9| . Assuming a density for the condensate of the form 

1 ( x y z \ 
n(x,y,z) = —— n , (18) 

where X r and X z are scaling variational parameters and 
no is a spherically symmetric density distribution (for in- 
stance the usual Thomas- Fermi inverted parabola). The 
condensate anisotropy \ r /\ z can be found by minimiz- 
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FIG. 3. Stability phase diagram in the ip-l plane 
(I — Uz/Lu r ) obtained in the Thomas-Fermi limit (N — > oo) 
and with £dd = 4. Condensate are always stable in the vicin- 
ity of the magic angle ifM = 54.7°. 



ing the expectation value of the total energy H to t 
H ho + H int where 



H ho = ciNm (2w 2 A 2 + uo 2 z \ 2 z ) 
is the harmonic potential energy and 



C2 



N 2 g 
X 2 X Z 



' 3 cos 2 ip — 1 \ / A, 
e dd{ / - 



(19) 



(20) 



the mean-field interaction energy, where c\ and C2 are 
numerical constants [EOl and 



/(«) = 



1 + 2k 2 3K 2 tanh _1 VT 



1 



(1- K 2 ) 



3/2 



(21) 



derives from an angular integration. 

The stability can be varied for a fixed value of £dd and 
the trap anisotropy I just varying one easy experimental 
parameter, namely the angle p. In Fig. 3 we show one 
example of how the stability phase diagram for e d d = 4 
can be varied as function of the angle p. For very large 
£dd the stable region in the p-l plane becomes a thin 
region around the magic angle line. 

Conclusion — We have shown that the anisotropy of 
the dipole-dipole interaction of a polarized gas can be 
used to tune the strength of the dipolar coupling and to 
excite collective excitations. Tuning the dipolar coupling 
will become an important tool for designing atomic quan- 
tum gases with novel properties and dynamical aspects. 
This approach in optically trapped Bose condensed gases 
should allow to make even very weak dipolar coupling vis- 
ible with experimentally feasible parameters (like B fields 
in the Gauss region and spinning frequencies up to sev- 
eral 10 kHz). The tuning of the interaction might become 
a relevant technique for stabilizing strongly interacting 
dipolar systems like heteronuclear molecular gases. 
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